Superconducting PrOs4Sbi2: a thermal conductivity study 
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The superconducting state of the heavy fermion PrOs4Sbi2 is studied by heat transport measure- 
ments on a highly homogeneous single crystal exhibiting only one transition peak in the specific 
heat. The field and temperature dependence of the thermal conductivity confirm multiband super- 
conductivity and point to fully open gaps on the whole Fermi surface. 



PACS numbers: 71.27. +a, 74.25.Fy, 74.25. Op, 74.45. +c, 74.70.Tx 
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Several unusual features have been reported on the 
superconducting state of PrOs4Sbi 2 , the first Pr-based 
heavy fermion (HF) superconductor |l|: a double super- 
conducting transition in the specific heat has been ob- 
served on different samples 0, IE El , like in the well- 
known case of UPt3, but its intrinsic nature has not been 
clearly established yet 0, E|. Other experiments also 
point to unconventional superconductivity in this com- 
pound: thermal conductivity (re) measurements in a ro- 
tated magnetic field , London penetration depth stud- 
ies [U and flux-line lattice distortion support nodes 
of the gap. These results contrast with scanning tun- 
neling spectroscopy (STM) [Tof . Sb nuclear quadrupole 
resonance (NQR) [ll| or mu on spin relaxation (^iSR) 
which measured a fully opened gap. Our first very low 
temperature re measurements under magnetic field (sam- 
ple A) provided strong evidence for multiband supercon- 
ductivity (MBSC) in PrOs4Sbi2, but sample quality did 
not allow analysis of the gap topology from the low tem- 
perature regime of re 0. 

In this Letter, we report a new study of thermal trans- 
port at very low temperatures on a highly homogeneous 
PrOs4Sbi2 single crystal. Supplementary specific heat 
measurements of this sample show only one single, sharp 
superconducting jump at T c , supporting an extrinsic ori- 
gin for the double transition reported so far. Improved 
sample quality also has profound impact on thermal 
transport, mainly on the temperature dependence k(T) 
in zero field. It provides compelling evidence for a rather 
"conventional" MBSC scenario with fully opened gaps on 
the whole Fermi surface. 

Our thin, platelet-shaped PrOs4Sbi2 single crystal 
(sample B2, - 760 x 340 x 45 ^m 3 , T c ~ 1.75 K, resid- 
ual resistivity (p) ratio p(300 K)/p(T c ) ~ 30 instead of 
~ 15 in sample A) has been extracted (gently "grinded 
down" against the disk of a diamond saw) from a con- 
glomerate of several small cubes of PrOs4Sbi2 (sample 
Bl^lx0.75x0.6 mm 3 ), grown by the Sb-flux method 
0. Specific heat (C p ) in zero field has been measured 
on a PPMS, first the entire conglomerate (Bl) and then 
sample B2 alone (see fig. k(T,H) of sample B2 par- 
allel to the magnetic field has been measured in a di- 
lution refrigerator by a standard two-thermomcters-one 
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FIG. 1: Specific heat C p /T, electric resistivity p and thermal 
conductivity n(T)/T in zero field around T c , documenting the 
collapse of the double transition (observed only in sample Bl) 
and the high homogeneity of sample B2. Nevertheless, a kink 
at 1.9 K on p(T) reveals some "traces" of the upper transition 
in sample B2, but the bulk superconducting transition clearly 
corresponds to the lower T c in C P /T(T) of sample Bl. 



heater steady-state method down to 30 mK and up to 
3 T (p H c2 {T -> 0) ~ 2.2 T). The carbon thermometers 
were thermalized on the sample by gold wires, held by 
silver paint on gold stripes evaporated on the surface of 
the sample after ion gun etching. The gold stripes are 
essential for the stability and the quality of the electri- 
cal contacts (resistance R% « lOrnfi at 4 K). The same 
contacts and gold wires were used to measure the electric 
resistivity of the sample by a standard four-point lock-in 
technique. The reliability of the experimental setup was 
checked against the Wiedemann-Franz law, givin g qu an- 
titatively similar results for L/L as reported in |13| . 

The excellent homogeneity of our sample B2 is docu- 
mented by the fact that the bulk superconducting tran- 
sition appears at exactly the same temperature on C p , 
k/T and p (see fig. ^ and below). Another criterion 
regarding crystal purity is the residual value of k/T 
in the T — > limit. For platelet B2, it is smaller 
than 1.6/tWK~ cm -1 , which corresponds to 0.07% of 
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FIG. 2: k/T(T) in zero field: comparison of former (A, open 
squares) and present (B2, open circles) sample. For sample 
B2, we separate phonon K ph (dash-dotted line) and electronic 
K el (crosses) contribution. The dashed line is a fit of K el within 
the MBSC scenario, see text. The inset shows the estimated 
temperature dependence of the phonon mean free path. 



k/T(T -> 0,p a H = 2.5 T > H c2 ) and is significantly 
lower than in former sample A. These signatures of high 
sample quality allow us to use thermal transport at very 
low temperatures as a sensitive probe of the low lying 
energy excitations in PrOs4Sbi2- 

Figure ^ displays the superconducting transition at 
zero magnetic field as seen by specific heat for samples Bl 
and B2, and in addition for sample B2 the corresponding 
k/T and p(T) curves. The specific heat of sample BI ex- 
hibits a double transition, comparable to those reported 
in The fact that both transitions behave similarly 
under magnetic field, and that the upper one always ap- 
pears inhomogeneous Q had cast doubt on the intrinsic 
nature of this double transition. The remarkable result 
is that the double transition collapses to a single, sharp 
jump of C p (at the lowest T c and of about the same overall 
height) just by reducing the crystal dimensions (sample 
Bl — > B2). Obviously, areas with a single and a double 
transition coexist within the same sample. It suggests 
that, like in URu2Si2 Q, the observed double super- 
conducting transition in PrOs4Sbi2 is related to sample 
inhomogeneity. A hint for extracting samples with a sin- 
gle transition comes from the preparation stage: in order 
to remove all small cavities appearing during the sawing 
process of Bl, we had to reduce the thickness of B2 down 
to only 50 pui. Similarly, such tiny dimensions were re- 
orted for another sample exhibiting a single C p jump 
. Further systematic (structural) investigations are re- 
quired to determine the nature of defects which might be 
at the origin of the sharp double transition in PrOs4Sbi2. 
In figure[21 we compare the temperature dependence of 



k/T in zero field of samples A (former) and B2 (present). 
For B2, k/T has dropped by about two orders of mag- 
nitude from T c down to 30 mK, with a clean T 3 behav- 
ior of k below 100 mK. In sample A, the low tempera- 
ture behavior of k/T is probably dominated by impuri- 
ties, defects and/or other inhomogeneities, resulting in a 
sort of cross-over regime k/T ~ T 0|. By contrast, at 
higher temperatures, k/T is qualitatively the same for 
both samples, except that there is now a sharp peak in 
k/T precisely at the superconducting transition (as seen 
by C p and p). In sample A, these features were much 
broader and the local maximum of k/T appeared below 
the resistive T c or the onset of the specific heat jump. 

In figureini we compare the normalized k(H ) /T data at 
50 mK of samples A and B2: in this temperature region, 
the quasiparticle mean free path is governed by elastic 
impurity scattering 0. The very fast increase of ther- 
mal conductivity at low fields is perfectly reproducible, 
and even more pronounced in sample B2 because of the 
significant drop of k/T for T <C T c in zero field (see 
fig. EJ: a magnetic field of only 5% of £f C 2(0) is enough 
to restore about 40% of the normal state heat transport. 
This robust feature, similar to observations in MgB2 , 
confirms MBSC in PrOs 4 Sbi 2 HH]. The plateau at 0.4k„ 
can be interpreted as the "normal state" contribution of 
the small gap band at T — > 0, observed once the vortex 
cores of that band completely overlap. Moreover, despite 
improved sample quality, testified by the sharp change of 
slope at T c , there is clearly no sign of a phase transition in 
the mixed state as suggested in 0: k(H,T^ 0) has no 
anomaly at the proposed H* line (H*(T ->0)« 0.8 T), 
whereas the B — > C transition in UPt3 was clearly seen 
on k{H) 0). 

Because of improved sample homogeneity, it is now 
worth analyzing the temperature behavior of k/T. Just 
above T c , L/L < 1 0, which indicates a phonon ther- 
mal conductivity (k p ' 1 ) negligible compared to the elec- 
tronic heat transport (k c1 ) in the neighborhood of T c . 
The change of slope observed at T c (d(ns / kjv) / d(T /T c )) 
is of order 1.4. In conventional superconductors, it is 
generally ascribed to the combined effects of the opening 
of the gap and the energy dependence of the electron- 
phonon scattering rate on K ei [l7l lisj . In the BCS 
weak-coupling limit, its maximum value is of order 1.4, 
when lattice scattering is the limiting mechanism for 
K el (see measurements on very pure In 0, [2(1 El)). 
For PrOs4Sbi2, electronic inelastic scattering may re- 
place the effect of electron-lattice scattering. Neverthe- 
less, taking into account the relative weight of elastic 
to inelastic scattering, as well as MBSC (negligible ef- 
fect of gap opening in the small gap band), the value of 
d(ns / kjss) /d(T /T c ) as 1.4 appears very large. It is likely 
a signature of strong-coupling effects, as observed (and 
calculated) for example in lead (d(Ks / kjv) / d(T /T c ) as 7 
EUgl). Indeed, Sb NQR [HI or heat capacity anal- 
ysis |5J have already stressed strong-coupling effects in 
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FIG. 3: k/T(H) (normalized to its value at the superconduct- 
ing transition) at 50 mK for samples A and B. The arrows 
indicate the lower (-ffci) and upper (H C 2) critical field. The 
data in the "field cooled" mode reveal residual flux pinning 
below 50 mT and a sensitivity of k to fields as low as 5 mT. 



PrOs 4 Sbi 2 . 

As regards now the origin of the local maximum of 
k/T slightly below 1 K, it had been ascribed to enhanced 
inelastic mean free path due to the condensation of elec- 
tronic scattering centers for T < T c . But the question 
of phonon or electronic origin (or both) was left open 
|l3|. It is now possible to give a quantitative estimate 
of the phonon contribution: indeed, a maximum of k/T 
below T c followed by a T 3 behavior of k at low tem- 
perature is well-documented from superconducting Pb 
and Nb [U [2j| , from the rare earth nickel borocarbides 
RNi 2 B 2 C (R=Lu, Y) EE 113 and man y other mate- 
rials. 

On cooling the phonon mean free path l ph increases 
from a law l ph ~ T _1 (when it is limited by electron- 
phonon interactions above T c ) up to a typical crystal di- 
mension at the lowest temperature. But an intermedi- 
ate regime, described empirically by l ph ~ T~ 1 (T C /T) n : 
starts below T c due to the reduction of scattering by elec- 
trons (see inset of fig. |2J). Quantitatively, an estimation 
for K ph /T is plotted in fig. 0from: 
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formal = with a = 60 MW/K 3 cm (from de- 
viation of the Lorenz number above T c |l3j). blf is fixed 
by k/T at lowest temperatures (T < 100 mK), yielding 
if ~ 60 /im (b = 10.9 x 10 3 WKT 4 m- 2 ). if is of order 
the smallest crystal dimension. The adjustable parame- 
ter is mainly the power law (n) for the boosted tempera- 
ture dependence of l ph : it proved impossible (adjusting 
n) to account for the local maximum only by the phonon 
contribution, so that a situation similar to that of Pb [22j 



or the high-T c oxydes |2jj, |28|, |29j is recovered. In any 
case, it is seen that K ph should follow a T 3 behavior up 
to 0.3 K, giving, in the temperature range 0.1-0.3 K, a 
robust estimate of the electronic contribution K el /T. 

Nevertheless, we tried to understand the normalized 
electronic contribution K el /t| 5x up to T < 0.6 K (region 
with dominant elastic impurity scattering). The most 
striking feature is that it is not possible to fit K el (T) if 
one assumes a BCS gap corresponding to T c = 1.729 K: 
K el /T starts to rise at much lower temperatures than 
expected, requiring a smaller gap value. Of course, this 
cannot be compensated by strong coupling effects (which 
only make it worse, increasing the ratio A/T c ), nor by 
another estimation of the phonon contribution (k pH (T) 
cannot be larger than blfT 3 , which is constrained by 
the measurements below 0.1 K). The data points can 
be quantitatively reproduced within a MBSC scenario, 
when we include a small A S (T) and a large Aj(T) gap 
function with the same T c , and two associated conduction 
channels: K el /T = n s ■ K e ljT+ (1 -n s ) ■ k^ /T. The best 
data fit is then obtained for a zero temperature gap ratio 
of about Ai/A S (T -> 0) ~ 3 with A S (T -> 0) ~ IK, 
and a "weight" for the small gap band n s ~ 0.35. This 
value is close to the 40% deduced from the "plateau" of 
k(H,T -> 0) (fig. El). The characteristic field scale Bf 2 
for the vortex core overlap of the small band gap can 

now be estimated from H c i/H^ 2 ~ ( A ' Z*f ) , where vf 



is the average Fermi velocity of band i. In [4J we assumed 
that the small gap band is also a light carrier band, with 
vf,s/vf,i ~ 5. We then get H^ 2 ~ 10 mT, which is of 
the order of H c i and seems reasonable owing to the data 
of k(H). So the main outcome of the analysis of this 
new data is the existence of a small but finite gap A s 
in PrOs4Sbi2, quantitatively consistent with the MBSC 
scenario deduced from k(H). 

Various measurements on the superconducting state 
of PrOs4Sbi2 have been interpreted either as poi nting 
to gap nodes 0, H @ or fully open gaps [H El 
Focusing on the latest, we can compare the extracted 

gap values. The NQR [ill as wel1 as the A* SR El stud- 
ies propose large ratios of 2A/fc B T c , respectively ~ 5.2 
and ~ 4.2, supporting strong coupling effects but not 
the presence of a small gap. Nevertheless, the NQR data 
shows a large residual relaxation rate (l/7i) below 0.5 K, 
which may point, as for our sample A, to crystal inhomo- 
geneities which prevent observation of the smallest gap. 
Moreover, like specific heat, the nuclear relaxation rate 
should be rather sensitive to bands with large density of 
states. So if our interpretation E3 of the small gap band 
as being also a "light mass" band is correct, it may have 
indeed little contribution to l/T\. The muon relaxation 
rate (a s ) measured with /iSR is controlled by the field dis- 
tribution which may not put more weight on the heavy 
than on the light bands. But the measurements were per- 
formed in a residual field of 20 mT, already larger than 
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FIG. 4: Electric R%{T) (full symbols) and thermal 
Rl h (T).L T (open symbols) resistance of the Au-PrOs4Sbi 2 
or Au-CeCoIns contacts. Full lines: fits of Rc h (T) according 
to equation R e c (T) are temperature and field indepen- 
dent, but Rl h (T) shows a highly singular behavior in zero 
field below 200mK on PrOs4Sbi2. This further supports a 
fully open gap in this system (see text). 



H^ 2 : so that again a s (T) is probably governed by the 
high energy excitations. But the "unusual" non linear 
field dependence of a s (T = 0.1 K) compares well with 
k(H): the MBSC scenario, with the small gap band as- 
sociated to light carriers can even "explain" the increase 
of a s at low fields as a s oc 1/m*. Eventually, the STM 
measurements ^3 proposed a gap distribution, which 
may extend from 120 /zV to 325 fiY (2A/k B T c -1.5-4.1), 
not so far from our analysis of n(T) (2As /fcsT c — 1.15, 
2A L /k B T c -3.5). 

A further robust and reproducible |l3| experimental 
observation supports a fully open gap in PrOs4Sbi2. In- 
deed, we can also measure on our setup the electrical con- 
tact resistance {R e c {T)) between the sample and the gold 
wire of the thermometer thermalisation, as well as its 
thermal contact resistance il*" is defined as P/AT, 
with P the Joule power dissipated in R e c {T) and AT the 
thermal gradient produced by P. For large contact areas, 
we expect that at low temperature 
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(aT 2 being the phonon contribution to the thermal con- 
ductivity of the contact). Figure0]shows that this is well 
observed below 1 K for PrOs4Sbi2 in a field of 20 mT, or 
CeCoIns in zero field (new measurements on this same 
setup). But for PrOs4Sbi2 in zero field, an unexpected 
divergence is observed below 200 mK, with no correspon- 
dence on R%{T) (field and temperature independent). On 
the same footing, CeCoIn5 in zero field, like PrOs4Sbi2 in 
20 mT, could be cooled below 10 mK whereas PrOs4Sbi2 



in zero field remained stuck above 25 mK. A reason 
could be that electronic heat transport is suppressed at 
the normal-superconducting interface when k B T <C A s 
(whereas electric current can always go through thanks 
to Andreev processes 01 course, this barrier is sup- 

pressed in very low field in PrOs4Sbi 2 (multiband effects) 
or in zero field in CeCoIns which has line nodes of the 
gap 0. 

In conclusion, we measured the thermal conductivity 
of a highly homogeneous PrOs4Sbi2 single crystal ex- 
hibiting a single jump on C p at T c . The reproducible 
field dependence k(H) at T < T c confirms the proposed 
MBSC scenario. Further support now comes from the low 
temperature k(T) and i?*' l (T) data which both point to 
isotropic, fully opened gap functions, in comparison with 
other measurements (NQR, uSR, STM). Owing to the 
still mysterious homogeneity problems in this compound 
and to the strong field sensitivity, analysis of the data re- 
quires a close look at the experimental conditions (sam- 
ple, field and temperature range) which may explain the 
remaining divergences of interpretations. 
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